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We present results for the electronic structure of plutonium by using a recently developed 
quasiparticle self-consistent GW method (QSGW). We consider a paramagnetic solution with- 
out spin-orbit interaction as a function of volume for the face-centered cubic (fee) unit cell. 
We span unit-cell volumes ranging from 10% greater than the equilibrium volume of the S 
phase to 90 % of the equivalent for the a phase of Pu. The self-consistent GW quasiparticle 
'"^ ■ energies are compared to those obtained within the Local Density Approximation (LDA). 

pH ■ The goal of the calculations is to understand systematic trends in the effects of electronic 

correlations on the quasiparticle energy bands of Pu as a function of the localization of the 
/ orbitals. We show that correlation effects narrow the / bands in two significantly different 
ways. Besides the expected narrowing of individual / bands (flatter dispersion), we find that 
an even more significant effect on the / bands is a decrease in the crystal-field splitting of the 
different bands. 

A^y^ • Keywords: Plutonium, electron correlations, GW approximation, first-principles electronic 

' structure 

oo 

' 1. Introduction 

o\ ■ 

' Much of our modern understanding of electronic correlations in narrow-band sys- 

tems has derived from many-body treatments of model-Hamiltonian systems such 
as the Hubbard model and the Anderson model. For example, with respect to plu- 
^ I tonium in particular, dynamical mean-field theory (DMFT) approaches have been 

very useful in elucidating the physics of the very strong correlations in this material 
(see, for example, Refs. [V\ and [2\ and references therein). Nonetheless, these calcu- 
lations are not first principles, and most of the physics comes from the model part 
of the Hamiltonian rather than the band-structure part of the calculations. Thus, it 
is still important to better understand the multi-orbital and hybridization effects in 
more realistic electronic-structure approaches that are less model dependent. The 
GW method is our best modern tool to examine these effects, because it includes 
correlation effects beyond that of conventional local-density approximation (LDA) 
band-structure techniques and yet is still first-principles. 

In this paper we study correlation effects of fee Pu as a function of volume. Our 
goal is not to specifically elucidate the correlation physics of Pu itself, since the GW 
approach is a low-order approximation and cannot treat the very strong correlations 
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of the 5 phase of Pu. Rather, we wish to use the volume dependence to tune the ma- 
terial from high atomic density (high pressure), where correlations can be greatly 
reduced due to the large hybridization between the Pu / orbitals, to low atomic 
densities (where the pressure would actually be in tension), where the correlations 
effects are very strong. By following this procedure, we can understand how corre- 
lations modify the properties of a material within the GW approximation. When 
the correlations become strong, it is certainly the case that higher-order approxi- 
mations like DMFT are necessary to accurately describe the material. Nonetheless, 
since GW is probably the correct starting point for such types of more sophisticated 
approaches, it is still useful to understand what happens to the electronic-structure 
of the material within the GW approximation as a function of the strength of the 
correlation effects, and, in particular, their effects on shifts in quasi-particle ener- 
gies, which is what GW is best at representing. 

In order to focus more specifically on the effects of correlations, we ignore one 
significant aspect of the electronic structure, viz., the spin-orbit coupling, even 
though this is important for a detailed comparison with experiment. Spin-orbit 
coupling mainly shifts / states around in energy, but at the same time add to the 
complexity of the individual / bands, hence tending to hide some of the correlation 
physics in a bewildering array of bands. Spin-orbit effects can easily be added in 
when a more accurate comparison with experiment is desired (as was done for an 
earlier paper on less correlated a uranium [3|])- 

For the same reasons we also ignore well-known large changes in crystal structure 
of Pu metal with volume, and present calculations only for the simple fee crystal 
structure as a function of volume. We study a range of Pu atomic volumes extending 
from well below that pertinent to the ground state a phase to well above that of 
the high-temperature 5 phase. The actual crystal structure of aPu is a complicated 
monoclinic structure with 16 atoms per unit cell, while 5Pu has the fee structure. 
Since JPu is well known to be a strongly correlated-electron metal while aPu 
appears to be reasonably well treated by conventional band-structure methods, 
we believe that our range of volumes corresponds to tuning the correlation effects 
between weak to moderate (small volumes) to strongly correlated (large volumes). 

In the calculations to be presented in the following we mainly focus on changes 
in the effective bandwidth of the / states in Pu. We will show that crystal-field 
effects are actually a more important factor in determining this bandwidth than 
the expected change in dispersion (flattening of the bands). 

2. Method 

The GW approximation can be viewed as the first term in the expansion of the non- 
local energy-dependent self-energy E(r, r',a;) in the screened Coulomb interaction 
W . Prom a more physical point of view it can be interpreted as a dynamically 
screened Hartree-Fock approximation plus a Coulomb-hole contribution [3]. It is 
also a prescription for mapping the non-interacting Green function onto the dressed 
Green's function: G^ — > G. This prescription can be described as follows. From the 
Hamiltonian 



(we use atomic Rydberg units: h = 2m = = 1, where m and e are the 

mass and charge of the electron) G^ = {uj — Hq)~^ may be constructed. Often 
G^ is calculated from the LDA eigenvalues and eigenfunctions; however, there is 
no formal restriction for how to choose the initial starting point G^. Then, using 




(1) 
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the Random Phase Approximation (RPA) 0], we can construct the polarization 
function D and screened Coulomb interaction W as 

D = -iG° X G° (2) 

and 

where u is the bare Coulomb interaction. The new Green's function is defined as 

^ " w - (- + V'^-' + + ' 
where V^^^ is the potential due to ions (Madelung) and is the Hartree potential 

FH(r) = 2 [ dr'^^ . (5) 

The single-particle density n(r) is calculated as n(r) ~ (ici;G'^(r, r, Li;)e*'^'^. 

In addition to this mapping of —>■ G, one can also generate an excellent 
effective potential from G that makes it possible to approximately do the inverse 
mapping of G ^ GO 0. The qSGW is a method to specify this (nearly) optimal 
mapping of G — > G", so that G" — > G — > G^ — > ... can be iterated to self-consistency 
0, [^. At self-consistency the quasiparticle energies of G^ coincide with those of 
G. Thus QSGTy is a self-consistent perturbation theory, where the self-consistency 
condition is constructed to minimize the size of the perturbation. The QSGW 
method is parameter- free, and independent of basis set as well as the LDA starting 
point & 7[. We have previously shown that QSGW reliably describes a wide range 
of spd [5, §-ll|, and rare-earth systems. We have also applied the method to 
calculate the electronic structure of a- uranium [§] . 

Our version of the QSGW method is based upon the Full Potential Linear Muffin 
Tin Orbital (FP-LMTO) method [IS], which makes no approximation for the shape 
of the crystal potential. The smoothed LMTO basis ||7] includes orbitals with / < 
^max = 6; both 7p and 6p as well as both 5/ and 6/ are included in the basis. The 
6/ orbitals are added as a local orbital |7j], which is confined to the augmentation 
sphere and has no envelope function. The 7p orbital is added as a kind of extended 
'local orbital,' the 'head' of which is evaluated at an energy far above Fermi level 
[3] and instead of making the orbital vanish at the augmentation radius a smooth 
Hankel 'tail' is attached to the orbital. The 7p and 6/ orbitals are necessary to 
obtain an accurate description of highlying bands, which are important for the 
accuracy of the polarization function in Eq. ([2]). For our calculations we use the fee 
Pu lattice with the following lattice constants a=4.11 A (which corresponds to 90% 
of the a-Pu equilibrium volume), 4.26 A (at the q-Pu equilibrium volume), 4.64 
A (at the J-Pu equilibrium volume) and 4.79 A (at 110% of the 6-Pu equilibrium 
volume) . 



3. Results 



In Fig.l we compare the QSGW one-particle electronic structure of 6-Fu with 
the LDA band-structure results. The spin-orbit interaction is not included in this 
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QSGW LDA 




Figure 1. (color online), (a) The QSGW energy bands (or quasi-particle energies) for (5-Pu along two 
symmetry directions (left panel), compared to (b) the LDA energy bands (right panel); the MuUiken 
weights of the / orbitals are presented in red (dark gray ) , for s orbitals in black and for d orbitals in blue 
(light gray), (c) Comparison of the total density of states (DOS) for QSGW, red (dark gray) solid line, 
and LDA, blue (light gray) dashed line. The Fermi energy is set at zero. 



calculation. The Mulliken weights of the / orbitals are presented in red (dark 
gray), of the s orbitals in black and of the d orbitals in blue (light gray). In both 
cases, the narrow bands located between -2 and 2 eV are predominantly due to 
seven 5/ orbitals. At the F point they are split by the cubic crystal field into one 
nondegenerate and two three- fold degenerate states. This degeneracy is reduced at 
general k-points in the Brillouin zone. The lowest dispersive band centered around 
-4 eV has primarily s character and the unoccupied bands above 2 eV are mainly 
due to d orbitals. At the T point the five d states are split by the cubic crystal 
field into one two-fold degenerate and one three-fold degenerate state. The sdf 
hybridization along the symmetry directions presented in Fig. [T] is generally very 
weak except for near the X point along the T — X direction, where there is very 
strong sd and df hybridization for some of the d and / branches. The degree of the 
hybridization is very similar in both the LDA and QSGW calculations, as is the 
center of all of the bands. The largest visible change is a significant narrowing of the 
5/-band complex. This effect has two components: first, the crystal- field splitting 
of the / bands is significantly reduced in QSGW, second, the bandwidth of each 
individual / branch (band) is also reduced in QSGW. The total effect appears as 
an overall narrowing of the total density of states (DOS) around the Fermi level 
(see Fig. (He)). In addition, since the area under the curve is proportional to the 
number of 5/ states, which remains constant, the amplitude of the quasiparticle 
peaks are also higher in QSGW. For example, the total DOS at the Fermi level in 
QSGW is 23 states/eV while in LDA it is 9 states/eV. 

The partial DOS is presented in Fig. [2j In both calculations the partial-5/ DOS is 
concentrated in a narrow energy interval around Ep. The partial-s DOS is mainly 
located in the occupied energy spectrum between -5 eV and Ep and the d bands 
are spread in a wide energy interval in both the unoccupied and occupied part of 
the spectrum with a few pronounced peaks at various energies. Overall, the QSGW 
and LDA s and d peaks are located at the same energies, with the exception of 
the narrow peaks around Ep, in which case the QSGW are visibly shifted closer 
to Ep. In this region these bands can be highly hybridized with the / states, and 
thus these shifts reflect the band narrowing of the / states. The bottom of the s 
and d bands relative to the position of / band is approximately at the same energy 
location in QSGVl^ and LDA. The / occupation changes in QSGW from its LDA 
value of 5.06 to 4.85 states. A similar reduction was observed in our calculations 
for a- uranium [31. 
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In Fig. [3] we present side by side the band structures for two different fee volumes. 
On the left side is the band structure and partial s and d DOS for the a volume 
(weakly correlated) and on the right side is the band structure and partial s and d 
DOS for the 6 volume (strongly correlated). In all cases, the red solid lines represent 
the QSGW results and the dashed-blue lines show the LDA results. Two major 
effects are seen in this plot: (1) the / bands narrow considerably with expanded 
volume, and (2) there is a similar effect on the d bands (notice the downward shift 
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Figure 2. (color online). Comparison of QSGW, red (dark gray) solid line, and LDA, blue (light gray) 
dashed line, partial DOS for (a) the d orbitals, (b) the / orbitals, and (c) the s orbitals. 
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Figure 3. (color online). Comparison of the QSGW, red (dark gray) solid line, and LDA, blue (light gray) 
dashed line, band structure. On the left side are the energy bands for the more weakly correlated case for 
a=4.26 A (equivalent to the density of atoms for the o-Pu equilibrium volume), and on the right side the 
more strongly correlated case of a=4.64 A (the 5-Pu equilibrium volume). 
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L r X 

Figure 4. (color online). Comparison of QSGW, red (dark gray) solid line, and LDA, blue (light gray) 
dashed line, energy bands along two symmetry directions. The magnitude of crystal field splittings Ai, A2 
are 1.02, 1.37 eV in LDA and 0.73, 0.57 eV in QSGW. A2 is reduced significantly in QSGW calculation. 
For the 5/ state at the F point, marked with a large dot, in Fig. |5]we present the energy dependence of 
the self-energy and calculated spectral function A{uj) 



of the bands at the highest energy as one goes from the a to the 5 volume). 

To examine band narrowing effects in more detail, in Fig. S] we expand the view 
of the 5/ bands. Also, in Fig. [5] we show the self-energy and spectral function 
for one of the QSGW5f states. Despite of the complicated energy dependence 
of the real and imaginary parts it appears that the state is described perfectly 
well by Landau's Fermi liquid theory. The imaginary part of the self-energy goes 
through zero at Ep and the spectral function has a single well defined peak centered 
at the QSGTy eigenvalue. There are no other pronounced features in the energy 
dependence of the spectral function. This is representative of all 5/ states. So 
the 5/ states within the QSGtytheory are well defined quasiparticles with very 
large lifetime around Ep; the QSGiyeigenvalues coincide with the quasiparticle 
energy. Therefore, our discussion is focused on the effects of electron correlations 
on the QSGiy eigenvalues. At the F point the / orbitals are split by the cubic 
crystal field into a nondegenerate state A2 and two three fold degenerate states 
Ti and T2. The splitting, Ai, between the nondegenerate state and the lowest of 
the three-fold degenerate states is 1.02 eV in LDA and 0.73 eV in QSGW. The 
splitting, A2, between the two three-fold degenerate states is equal to 1.37 eV in 
LDA but only 0.57 eV in QSGW. This significant reduction of the crystal field 
splitting in QSGW is the major part of the band narrowing observed in the DOS 
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Figure 5. The energy dependence of the real and imaginary part of self-energy together with the spectral 
function for the quasiparticle at k = (0, 0, 0) and i?o = —0.917 eV. 
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Table 1. The band width of Pu 5/ 
bands along L — T and P — X sym- 
metry directions. The band width 
is defined as the difference between 
the maximum and the minimum 
energy of a particular band along 



the direction. 


The width 


is given 


in cV 








L - r 


r - X 


A2 LDA 


0.48 


0.64 


A2 QSGW 


0.43 


1.21 


Ti LDA 


0.71 


0.83 


Ti QSGW 


0.38 


0.44 


T2 LDA 


0.77 


1.02 


T2 QSGW 


0.32 


0.45 



in Fig. [TJ Another important aspect is the reduction of the width of each individual 
/ band. In Table [1] we present the values of the band width of Pu 5/ bands along 
L — r and T — X symmetry directions. The band width is defined as the difference 
between the maximum and the minimum energy of a particular band along the 
direction. In all cases the band width is reduced significantly in QSGW. The A2 
band along F — X is a striking exception from this rule. As we show in Fig. [6] this 
band strongly hybridizes with the d states. At the F point it is 100% / but at 
the X point it is predominately d. All the other bands shown in Fig. H] remain 80 
to 100% / throughout the symmetry directions shown in the figure. This explains 
the anomalous change in the width of this band for this particular direction. From 
Fig. [3] it is evident that, while in QSGW the width of the / bands are significantly 
reduced, the width of d bands are practically the same as for LDA. The A2 band 
at the F point has mainly / character and therefore moves upward from its LDA 
position due to the significant reduction of the / crystal field in QSGW, but at the 
X point mainly has d character and therefore remains approximately at its LDA 
position (only slightly lower due to the slight downward shift of the center of d 
band in QSGW). The cumulative effect is that in QSGW this band is stretched. 

In Tables [2] and [3] we show the values for the crystal field splitting of 5/ bands and 
/ band widths along L — T for several volumes of the unit cell. The band width and 
crystal field splitting of all bands is reduced as we move from the lower to higher 
volume case. This is a result of the reduction of the / band relative extend in the 
crystal. It is also seen that the QSGW crystal field splittings and band widths for / 
orbitals are always smaller than in LDA. We have also considered the band width 




k (2jt/a) 

Figure 6. (color online). QSGW, red (dark gray) solid line, and LDA, blue (light gray) solid line /-orbital 
weight of the A2 band along V — X direction. Also shown is the QSGW, red (dark gray) dashed line, and 
LDA, blue (light gray) dashed line d-orbital weight of the A2 band along the same direction 
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Table 2. The crystal field splitting Ai and A2 
of Pu 5f states at the P point. The QSGW val- 
ues for the splitting arc significantly smaller than 
those of LDA calculation. The energy splitting is 
given in cV 





0.9 Va 




Vs 


1.1 Vs 


^LDA 

^hsGW 
^1 


2.14 
2.09 


1.75 
1.67 


1.02 
0.73 


0.84 
0.51 


^LDA 
"2 


2.86 
1.96 


2.33 
1.36 


1.37 
0.57 


1.11 
0.41 



Table 3. The band width of Pu 5/ bands along L — 
P symmetry directions for different volumes of the unit 
cell. The band width is defined as the difference between 
the maximum and the minimum energy of a particular 
band along the direction. We also provide the full width 
at half maximum (FWHM) of the broadened /-partial 
DOS shown in Fig.[T] The width is given in cV 





0.9 Vo. 




Vs 


1.1 Vs 


A2 LDA 
A2 QSGW 


1.42 
1.23 


1.0 
0.94 


0.48 
0.43 


0.37 
0.30 


Ti LDA 
Ti QSGW 


1.38 
1.17 


1.15 
0.86 


0.71 
0.38 


0.6 
0.28 


T2 LDA 
T2 QSGW 


1.62 
1.30 


1.29 
0.84 


0.77 
0.32 


0.63 
0.23 


FWHM LDA 
FWHM QSGW 


1.28 
0.79 


1.09 
0.63 


0.78 
0.39 


0.72 
0.34 



of the entire 5/ band complex. To do this we have appUed a Gaussian broadening 
on the partial-/ DOS (Fig. [7]). In this case the 5/ band complex appears like a 
single large peak. We can define the width of this band as the full width at half 
maximum (FWHM) of the peak. This is also presented in Table [3l It is evident 
that in both calculations the width is reduced as the volume increases. The rate 
of reduction is the same in both calculations. But the QSGW FWHM is always 
significantly smaller than the LDA. Therefore, we conclude that in QSGW the / 
orbitals contract due to the more accurate treatment of correlations. On the other 
hand, the s and d bands are very itinerant and are therefore already described 
accurately at the level of the LDA treatment of correlations. 

4. Conclusions 

In conclusion, we have applied QSGW theory to 6 plutonium as a function of vol- 
ume in order to systematically understand the effects of electronic correlation on 
the band narrowing of the / bands. Unlike conventional model-Hamiltonian treat- 
ments of strongly correlated systems our approach is first principles and indepen- 
dent of any choice of model parameters, and hence provides a unique opportunity 
to examine the effect of electron correlations on the quasiparticle band structure as 
a function of /-orbital localization. In this way we have demonstrated that QSGW 
and LDA prediction for the s and d electron subsystems are quite similar. This 
is because these electrons are very itinerant and therefore their description lies 
within the validity of LDA. However, the QSGW 5/ bands are much narrower 
than their LDA counterpart. We believe that our results for the first time show 
significant details of the band narrowing due to electron correlations that have 
not been previously studied. In particular, the QSGW calculations show that the 



April 30, 2009 17:8 Philosophical Magazine pu-v5 



Philosophical Magazine 



O 
Q 



broadened 



r Vg+O.lVg ll 


. ' 1 

— A 


1 

- 1 


' 1 ' 
" 1 


1 1 1 


: 1 -t-r 




" V 

a 




'■ — 1 '1 "1 

: ' 1 ' i-n^ 





L V -O.IV 

a a 









E-Ep (eV) 



E-Ep (eV) 



2 -1.5 -1 -0.5 0.5 1 1.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 



Figure 7. (color online). The QSGW and LDA /-partial DOS for four different volumes. The panels on 
the right side show the corresponding DOS on the left panel broadened with a Gauss function. 

major contribution to band narrowing is actually the reduction of the crystal- field 
splitting of 5/ states as compared to effects from a reduction in dispersion (flat- 
tening of the individual bands). This is a significant change in the character of 
the 5/ states and suggests the importance of using GW approaches as input to 
more sophisticated correlation approaches such as dynamical mean-field theories 
(DMFT). 
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